GLOBAL WEAK SOLUTIONS FOR SOME OLDROYD MODELS 



OLFA BJAOUI AND MOHAMED MAJDOUB 



Abstract. We investigate an evolutive system of non-linear partial differential 
equations derived from Oldroyd models on Non-Newtonian flows. We prove global 
existence of weak solutions, in the case of a smooth bounded domain, for general 
initial data. The results hold true for the periodic case. 



. 1. Introduction 



Contents 



2. Notations and Functional spaces 

3. Technical lemmas and statement of the main results 

4. Proofs of the main results 

4.1. Galerkin approximation 

4.2. Uniform estimates 

4.3. Strong convergence of v in L 2 

4.4. Equi-integrability of t n in L 2 

4.5. Strong convergence of t n in L 2 

4.6. Weak limits of non linear terms involving D(v N ) 
References 



1. Introduction 



1 
] 

5 
9 
9 
11 
13 

U 

16 
21 
24 



The Navier-Stokes equation is a model that describes the evolution of incompress- 
ible fluids with constant viscosity called Newtonian fluids. However, many experi- 
ments show that the viscosity of a fluid may vary with the pressure, see Andrade 
[2] , see also the book by Bridgman [5] . Further details and references to more recent 
experimental studies can be found in the book by Szeri [35] and in the paper by 
Malek and Rajagopal [21]. There are also many experiments which show that the 
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viscosity may depend on the symmetric part of the velocity gradient. We can refer 
to Schowalter [31], Huilgol [IB] . Recently, Malek, Necas and Rajagopal [21], Bulicek, 
Majdoub and Malek [BJ have established existence results concerning the flows of flu- 
ids called non-Newtonian fluids, whose viscosity depends on both the pressure and 
the symmetric part of the velocity gradient. The constitutive equation for such fluids 
is given by 

(1.1) T = -pI + v(p,\D(v)\ 2 )D(v), 

where T is the Cauchy stress, p is the pressure, v is the velocity, D(v) is the symmetric 
part of the velocity gradient and v is the viscosity function. Thus, the evolution of 
such fluids is governed by the equation 

p(d t v + v.Vv ) = divT, 

where p is the density of the fluid. 

However, there are some fluids that do not obey to the constitutive equation (jl.ip 
such as Blood. Yet, it was shown experimentally that Blood is a complex rheological 
mixture that exhibits shear thinning and elastic behavior, see Thurston J3BJ- The 
constitutive equation for such fluids is given by 

T=-pI+ S, 

where S is the extra-stress tensor which is related to the kinematic variables through 
S + Ai^ = 2p(\D(v)f)D(v) + 2A 2 ^, 

H is the viscosity function, Xi > and A 2 > are viscoelastic constants. 
The symbol denotes the objective derivative of Oldroyd type defined by 

nc 

_ = d t S + v.VS + S.w(v) - w(v)S, 

where w(v) is the antisymmetric part of the velocity gradient. The extra-stress tensor 
S is decomposed into the sum of its Newtonian part 

r s = 2^D(v) 

and its viscoelastic part r e . The constitutive equation for r e is given by 

Dt p _ / „„, A? 



r e + A r 



2(p(\D(v)f)-^)D(v). 



Dt 

Hence, we get the generalized Oldroyd-B model given by 
pd t v + pv.Vv — ^Af + Vp = div(T e ), 

d t r e + J-r e + v.Vr e + r e .w{v) - w(v)r e = £ (Ml^)| 2 ) - %)d(v). 

Particularly, when p is constant, we recover the Oldroyd-B fluid with constant viscos- 
ity, see [28] . Existence of local strong solutions to the Oldroyd-B model was proved 
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by Guillope and Saut in jH] and [15]. Fernandez-Cara, Guillen and Ortega [TT], [T2] 
and [13] proved local well posedness in Sobolev spaces. In the frame of critical Besov 
spaces, Chemin and Masmoudi had proved local and global well-posedness results in 
[7J. In addition, non-blow up criteria for Oldroyd-B model were given in [T8] as well 
as in [21] . For the sake of completeness, we refer the reader to the most recent papers 
about Oldroyd models as [22 El E21 EH E5] . 

Global existence for small data was proved in [23] and [26]. Considering general 
initial data, Lions and Masmoudi [28J had established results of global existence of 
weak solutions when the viscosity function fi is constant. Our aim in this paper is to 
generalize the results in [28]. More precisely, we will focus on the following system 

' ' d t v + v.Vv - div (f(D(v))) + Vp = div r, in R+ x Q, 
d t r + w.Vr + t.w(v) — w(v).t + or — g(D(v)), in R + x fi, 
^ M div v = 0, in R+ x fl, 

k t»(0, x) = v (x), r(0, x) = r (x), 

where Q can be considered either the torus T n , or a smooth bounded domain of R n , 
n = 2, 3 and in this case (S) is supplemented by the Dirichlet homogeneous boundary 
condition. 

The function g(D(v)) is given by 

g(D(v)) = fl(\D(v)\ 2 )D(v) = ^(Ml - 0)D(v), 

where fi(\D(v)\ 2 ) = 1 - A + A(l + \D{v)\ 2 ) r -^ , X G [0,1], r e [1,2]. 

Physically, the parameters a and b are given by a = b = ^rp~, where We is the 

Weissenberg number and 9 e]0, 1[ is the ratio between the so called relaxation and 

times. 

Let us mention that the expression of fl is motivated by the system studied by Arada 
and Sequeira [3] in the steady case with / being the identity map, where existence 
of a unique solution was established for small and suitably regular data. Notice that 
when the term f(D(v)) is replaced by uAv, with v = -jj^, Re being the the Reynolds 
number, and r is equal to 2, then (S) turns into the system studied by Lions and 
Masmoudi in [28]. Let us remark that when b = and tq = 0, then r = solves the 
second equation in (S) which consequently will be reduced to the system studied in 
[SJ and [9]. 

Our objective is to prove global existence of weak solutions, for general initial data, 
under suitable hypotheses on the function /. 

The layout of this paper is as follows. In the following section, we give some notations 
and we introduce the functional spaces used along this paper. The third section is 
devoted to some technical lemmas and to the statement of the main results. In the 
fourth section, we prove existence of approximate solutions (V N ,T N ) to the system 
(S) by using Galerkin method. Then, we prove the equi-integrability in the Lebesgue 
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space L 2 ([0,T] x Q) of the sequence t n , and we derive an evolution equation for 
r] := \t n — t\ 2 . This enables us to get the strong convergence of the sequence t n in 
L 2 ([0, T]xfl). Therefore, we can identify the weak limits in the sense of distributions 
of non linear terms. 

2. Notations and Functional spaces 

Let us introduce some notations which we use throughout this paper. With M nxn 
we denote the real vector space of n x n matrices. We have (divr)i = £V djTij, for 
t G M nxn , and for a vector field v, we have (Vw)jj = djUi. 

For A = (dij) and B = (bij) in M nxn , A : 5 stands for the sum £V a^-ft^- and |A| 2 
denotes the quantity ^4 : A. 
Let us denote by: 

V = {<peV(Sl) n ,div<p = Q}, 

and 

V = {M e V(tt) nxn , M l = M}, 
when Q is a bounded domain of R", 

V = {0 g C^ r (fi) n , div<f> = 0, / 0cfe = 0}, 

and 

V = {M G C^ r (fi) nxn , M* = M, ( Mdx = 0} 

when Q is the torus T n , 

H = the closure of V in the L 2 (fi) n -norm, 

H = the closure of V in the L 2 (f2) raxra -norm, 

V 9 = closure of V in the L 9 (f2) n2 -norm of gradients, q>l. 

V g = closure of V in the L' 3 (f2) nxn -norm of gradients, q>l. 

V s = the closure of V with respect to the W"' (f2) n -norm, with s > 1 + |, 

Vs = the closure of V with respect to the VF"' 2 (f2) nxn -norm, with s > 1 + |. 

The condition on s is due to the fact that: if u G W s ' 2 (tt) n , then Vv G VF s_1 ' 2 (fi) n2 

and W 8 ' 1 ' 2 ^) L°°(fi) if ± - ^ < 0. 

We denote by V,' the dual space of V s and by <, > s the duality between V s and V!.. 
The scalar product in L 2 (f2) will be denoted by (.,.). 
For s > 1 + | and g > 2, we have the following inclusions 

Kcv,c^^'cv;c v?, 

and 

v;cv,ch-h'cv',c v' s . 

For a sequence /jv in V'([0,T] x fi), we denote by /at its weak limit in the sense of 
distributions. 
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By C we denote any constant that may depend on the n-dimensional Lebesgue 

measure of fl, v , r and on T > 0, but not on N. 

Let us mention that subsequences will not be relabeled. 

3. Technical lemmas and statement of the main results 

We give some technical lemmas needed for the proofs of our main results. 
One of the properties of the norm on Banach spaces is the lower semi- continuity given 
by the following lemma. For more details, we refer to [58] . 

Lemma 3.1. Let X be a Banach space equipped with a norm ||.||. Let (x^) be a 
sequence in X that converges weakly to some x in X , then (xjv) is bounded in X and 
we have 

llxll < liminf ||xjv||> 

iV-H-oo 

Lemma 3.2. (Vitali's Lemma) Let Q be a bounded domain in R ra and : — ?• R 

be integrable for every JVeN. Assume that 

(i) linuv^+oo f N (y) exists and is finite for almost all y G ; 

(ii) for every e > 0, there exists 5 > such that 

sup / \f N (y)\dy<e \/HcQ,\H\<5. 

Then 

lim / f N (y)dy= / lim f N (y)dy. 

For the proof of lemma I3.2[ we can refer to pQ . 

Lemma 3.3. (Aubin-Lions lemma) Let 1 < a,f3 < +oo and T > 0. Let X be a 
Banach space, and let Xo,Xi be separable and reflexive Banach spaces such that Xq 
is compactly embedded into X which is continuously embedded into X%, then 

{v E L a ([0,T),X o );d t v e L^([0, T},Xi)} is compactly embedded into L a ([0,T],X). 

For the proof of lemma l3T3| we refer to [27]. A generalized form of this lemma for 
locally convex spaces and j3 = 1 can be found in [33J. 

The following lemma plays an important role in the theory of existence of solutions 
to ordinary differential equations. The proof of such lemma can be found in 



Lemma 3.4. (Caratheodory Theorem) Let c : Is = [to — 5,t + 5] — > R™ the system 
of ordinary differential equations 

d c(t) = F(t,c(t)), tei s 
= c e R n , 

Assume F : 1$ x K — )■ R n ; where K = {c G R n , |c — c | < /3}, for some (5 > 0. 
If F satisfies the Caratheodory conditions: 
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(%) t (->■ Fi(t,c) is measurable for all i = 1 , n and /or all c E K, 

(ii) c4 c) is continuous for almost all t £ 7,5, 

fraj i/iere exists an integrable function G : 7,5 — >• R suc/i i/iai 

|Fi(t,c)| < G(t), V(t,c)eI 5 xK, Vz = l,...,n, 

t/ien t/iere exists 5' £]0, 5[ and a continuous function c : I& —> R n snc/i 

(i) exists for almost all t £ 

(ii) c solves E. 

The following theorem, the proof of which can be found in [I] , deals with compact 
injections of Sobolev spaces into Lebesgue ones. 

Theorem 3.1. (Rellich-Kondrachov) Let fl be a C 1 bounded domain ofH n . We have 
the following compact injections 

(i) ifp = n, then W lj, {Q) C L q (0), Wq £ [l,+oo[, 

(ii) ifp > n, then W 1 '^) C C(O). 

In order to prove the existence of approximate solutions, we use Galerkin method 
by considering a special basis in the space V S (Q) as in [8]. The proof of existence of 
such basis, see [30] , relies on solving the following spectral problem: find w r £ V s and 
A r £ R satisfying 

<w r ,ij> s =X r (w r ,ij), WipeV s . (V) 

Theorem 3.2. There exists a countable set {\ r }'^L l and a corresponding family of 
eigenvectors {w r }'^ =1 solving the problem (V) such that 

(i) (w r ,w r ')=5 r y Vr,r' £ N ; 

(ii) 1 < A x < A 2 < ... and X r — > oo as r — > oo, 

(m) <^,^=> s =5 ry Vr,r'£N ; 
(iv) {w r }'^ =1 forms a basis in V s . 

Moreover defining % := span{w l , ...,w N } (a linear hull) and 

N 

p N (v) ■= J2(v,w r )w r en N ,ve Vs, 
i=i 

we obtain 

II P II C(Vs,Vs) - Ij \\P \\c(V^,Vi) < 1, || -P \\c{H,H) < 1- 

The same results hold if we replace V s by V s . 
One of the tools that we use in order to prove our main results is Young measures. 
Thus, we should give some well-known facts about this tool. We can refer to [32] as 
well as to [17] and |30j for more details. 
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Theorem 3.3. Let (f N : Q — > R n ) be a bounded sequence in L 1 ^) and fx be its 
weak limit in V'(Q). Assume that (f^) is equi-integrable in then 



fN — ^dv x (\), a.ex G fi, 



where v is the Young measure generated by the sequence /n- 

Now, assume that fx converges weakly in L 1 (Q), then /jv is equi-integrable in L 1 (J7) 
and 



Jn= ^dv x (\), a.exEil. 



Theorem 3.4. Let (f^ : Q — > R n ) be a sequence of maps that generates the Young 
measure v. Let F : Q x R n H be a Caratheodory function that is a function that 
satisfies the Caratheodory conditions. Assume that the negative part F~(., f N (.)) is 
weakly relatively compact in L 1 ^), then 

/ / F(x, X)du x (\)dx < liminf / F(x, f N (x))dx. 

Now, we are ready to state our main results. The first Theorem gives a global 
existence result for the following system 



(Si) 



' d t v + v.Wv - div (f(D(v))) + Vp = div r, in R + x Q, 
d t r + jj.Vr + ar = g{D(v)), in R+ x Q, 
div v = 0, in R + x f2, 

y(0,x) = v (x), t(0,x) = t (x). 



Theorem 3.5. Let f : M nxn — > M nxn be a continuous function satisfying the fol- 
lowing hypotheses for some p G]2, +oo[ when n = 2 and p g]§, +oo[ when n = 3: 
(Hi) growth: there exists c > and c > such that 

\f(A)\ < 5 + c\A\ p -\ VA G M nxn , /(0) = 0, 

(H 2 ) monotonicity: there exists v > such that 

(f(A) - /(B)) : (A — B) > v{\A - B\ 2 + \A — B\ p ), WA, B G M nxn . 

Let vq and To be in "H(f2) and "H(f2) respectively. 

i) If v satisfies 2z/(l — 6) > 1, then for an arbitrary A G [0, 1], there exists a global 
weak solution (v,r) to the system (S\) such that 

v g L°°(R + ,n(n)) n l p (o,t, v p (n)), r g L°°(R + ,n(n)), vt > o. 

ii) If v is such that < 2z/(l — 9) < 1, then the same result holds for an arbitrary 
A G [0, y/2u(l-6)[. 



OLFA BJAOUI AND MOHAMED MAJDOUB 



Let us assume the existence of a scalar function U G C 2 (R r ' 
' /, such that for 
ij,k,l G {1, ...,n} 



of /, such that for some p > 1, C\ > 0, Gi > we have for all rj, £ G Af^/ 



called potential 

nxn and 



(3.2) 
(3.3) 
(3.4) 
(3.5) 



917(77) 
driij 



U(0) 



dU(0) 
drjij 



d 2 U(r]) 
dr/ijdrjki 



^H>Ci(i + Hr a |e| a , 



9 2 f/(77) 



driij dr]M 



<c 2 (i + \ v \y- 2 . 



As consequence of these assumptions, there exists C > such that 

\f(A)\<C(l + \A\y-\ 
Moreover, for p > 2, there exists ^ > such that 

{f{A) - f(B)) :{A-B)> v(\A - B\ 2 + \A- B\ p ), VA, B G M nxn . 
Standard examples of functions / whose potentials satisfy these assumptions are 

f(A) = (l + \A\r~ 2 A 

and 

f(A) = (l + \A\ 2 )^A. 

For more details about the existence of a such potential U and consequences of prop- 
erties E2J|33] we refer to [30] . 

The second Theorem concerns the following system 

' ' dtv + v.Vv - div (f(D(v))) + Vp = div r, in R+ x ft, 
d t r + v.Vt + ar = bD(v), in R+ x fl, 
div v — 0, in R + x £7, 

v(0,x) = v (x), t(0,x) = t (x). 

Theorem 3.6. Let f : M nxn — > M nxn be a C 1 -function satisfying the following 
hypotheses for some p G]2, +00 [ when n = 2 and p G]|, +00 [ u>/ien 72, = 3: 
(Hi) growth: there exists c > sitc/i that 

\f(A)\ < c\A\ p -\ VAeM nxn , 

(H2) coercivity: there exists v > such that 

f(A) : A > is\A\ p , VA G M nxn , 
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(H 3 ) monotonicity 

{f {A) - f(B)) : (A - B) > 0, VA, B e M nxn . 

Let Vq and To be in T-L{Q) and respectively. Then, there exists a global weak 

solution (v, t) to the system (S2) such that 

v G L M (fi + ,H(fi))nL p (0,T, V P (Q)), t G L°°(R + 1 "H(fi)), VT > 0. 

Obviously, if / satisfies the hypotheses of Theorem 13.51 with c = 0, then / satisfies 
the hypotheses of Theorem 13.61 

Let us remark that the quadratic term t.w(v n ) — w(v n ).t is not present in (S\) 
neither in (S2). The difficulty of this fact will be explained in the proofs of Theorems 
13.51 and 13.61 that will be given in the case of a bounded domain. However, they can 
be easily adapted to the periodic case. 



4. Proofs of the main results 

The results in subsections 14.11 14.2[ 14.31 and 14.41 hold for the system (S) and thus for 
the systems (S\) and (S2). 

4.1. Galerkin approximation. We will show the existence of approximate solution 
to the system (S) via Galerkin approximation as in 0. Hence, let {a r }^L 1 and {a r }^ =1 
be basis of V s and V s respectively given by Theorem 13.21 
Let T > and iV > 1 be fixed. We define 

N N 

v N (t,x) = Y,d^t)a k (x), r N (t,x) = Y,c%(t)a k (x), 

k=l k=l 

where the coefficients c^ (t) and d% (t) solve the so-called Galerkin system: a system 
of 2N nonlinear equations with 2N unknowns 

( j n d t v N <jdx + f Q (v N .V)v N o}dx + j n f{D{v N )) : Dtfdx = - J Q r N : Da'dx, 
(S N ) I J n 8 t r N : aHx + f n v n .Vt n : aHx + a J n r N : a j dx = J n fi(\D(v N )\ 2 )D(v N ) : a 
(cf (0) = (7-0,0^,^(0) = K,^),! < i,j<N 

The initial conditions on c% and are such that 

v N (0,x) = P N v (x), t n (0,x) = P n t (x), 

where P N and P N are the orthogonal continuous projectors of H and % respectively 
onto the linear hulls of the first eigenvectors a r , r = 1, ...N and a r ,r = 1, ...TV respec- 
tively. 
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The orthogonality of the two basis {a r } c ^° =1 and {a r }^L 1 in H and H respectively 
imply that the system (<SW) can be rewritten as 



±( r N ]N\ _ jr. .(f c n r N jN jN\ 

cf(0) = (r ,c^),di v (0) = Ka < ), 



where for 1 < i, j < N 

T (+ r' N r' N rl N J N \ — 

J ~i,j\ l -> c l > •••) c Afi tt l ' ■■•! a N ) — 



9. 

N 



J2 c k^ k ) : Da'dx- j (^rffa fe ).(^^Va fc )aVx 

k=l ^ Q k=l k=l 

„ N N 

f(J2^D(a k )) : D(a*)dx,- / dfa fc ).(^ cf Va fc ) : rfdz 
fe=i ^ n fc=i fe=i 

» W » N N 

a / (E c ^« fc ) : / mi^EWi^EW : a3dx ) 

J Q k=l J Q k=l k=l 

r N N 

^ Q fc=l fc=l 



A? 



A 



Let i? > and X C R 2JV be the ball of center (cf (0), 4(0), <(0), <(0)) and 
of radius R. We consider : [0,T] x K ->■ R 2JV . 

The continuity of /t and / lead to the continuity of J^j over [0,T] x X. In addition, 
thanks to the continuous inclusion W s ~ 1,2 (Vl) > L°°(Q) if \ — ^ < 0, we get 

. AT N 

^ n fc=i fe=i 
In the same way, we have 

N N N 

(Y,d^a k ).(Y, d k^a k )a^x\<C(R,N)\\a%2 {n) £ ||a fc || i2(n) .||Va fc '|| L oc (n) . 
fe=i fe=i fc,fe'=i 

From the growth hypothesis on /, we deduce that 

N , N 

fC£dkD{a k )):D{a*)dx\ < ~c\\D{a>)\\ Lm + c H^k D(a k )r l \D(a*) 
k=i ^ n u - 1 



fe=l 



AT _ 1 

< C(i?,iV)(^P(a fe )|| L . w ) P \\D(a*)\\ L - {a) 
k=i 



+ c||D(a*)|Uw 
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As fl is bounded, we estimate in the same way the remaining terms in T^j to get 

|^(t,cf,..., 4, <,...,<)| <C(R,N), 
where C(R, N) is a constant that does not depend on t. 

The standard Caratheodory theory provides the existence of continuous functions 
(cf , cjv, df , d$) solutions to (Sjv) at least for a short time interval with ^(c^, c$, 
is defined almost everywhere. The uniform estimates that we will derive in the next 
subsection enable us to extend the solution onto the whole time interval [0, T]. 

4.2. Uniform estimates. Multiplying the first equation in (S N ) by df(t), then tak- 
ing the sum over i — 1, N, we obtain 

Id, 

The monotonicity hypothesis on / leads to 

(4-6) —llt^Wlli^n) + v(\\D(v N (t))\\Un) + ll^(^(*)) 11^^)) < 

F r N : D(v N )dx. 



\v N (t)\\h (n) + [ f(D(v N )):D(v N )dx = - [ r N :D(v N )dx. 
Jn Jn 



J 

Jn 



o 

Multiplying the second equation in (SV) by Cj(t), then taking the sum over j = 
1, N, we obtain 

^411^)11^(0) + fl|r N (t)lli a( n) = r ^ / (Ml^)l 2 )-^:^). 
Notice that (t n .w(v n ) — w(v n ).t n ) : r N = since r N is symmetric. Hence, we get 

^11^(011^(0) + -\\D(v N (t))\\ P L ^ } + —ll^Wlli-cn) + ^II^W||! 2( o) 
< J^(\D(v N )\") - l)r N : D(v N )dx 



< 



Young inequality implies that 

(4.7^11^)11^^ + v\\D(v»(t))\\U«) + YbJt^ Nml ^) + f H^WII^W 
(4-8) < ^7^)^(^11^(0)11^) + ^11^)111,(0)). 
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Let 7 := 2 (i-fl) y^J - Notice that in both cases i) and n) in Theorem 13. h\ we have 
7 < 1. Notice also that to absorb the terms at the right-hand side of (4.8), we do 
not need that / satisfies the strong monotonicity condition. However, we need just 
that / satisfies the coercivity condition f(A) : A > u\A\ p . In fact, as ft is a bounded 
domain and p > 2, then thanks to Holder inequality and Young one 



\D(v" (t))\\Un) = / \D(v"(t))\'<b 
Jn 

< c{\n\)( /"(WO 



N, 

x Jn 

< C(\n\) + ^\\D(v N (t))\\ P LPi n) 

< c(\n\) + \\D(v N (t))\\ p Lm - 



Finally, we obtain 



(4-9) —H^WII^n) + «(l-l)\\D(v N (t))\\ p LP{n) + ~\\r N {t)\\h(p) 

a 



(4-10) + £(1 -7)11^)11^) <C 

Particularly, this leads to the fact that 

\(ci , c N , di , d^-)\ 2 < C, 

where C is a constant that does not depend on t neither on N. 

This uniform boundedness with the continuity of (cf\ c^, d±, d^), we deduce 

that the functions (c^ , ...,c^,di, —,d^) are defined on the whole interval [0,T], see 

Zeidler [39] for more details. 

By Korn inequality, we have 

\\Vv N (t)\\r» (m <C\\D(v N 



and we deduce that we can extract subsequences such that 

v N ^*v weakly in L°°(0,T,H) 
v N weakly in L p (0, T, V p ) 

t n ^* t weakly in L°° (0, T, H) 

t n ^ r weakly inL 2 (0,T,^). 

The growth condition on / and the uniform boundedness of the sequence D(v N ) in 
L p ([0,T] x ft) imply that f(D(v N )) is bounded in L p '([0,T] x ft), and thus 

f{D{v N )) f(D(v N )) weakly inZ/([0,T] x ft). 
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In addition, as g(D(v N )) is bounded in L P ([Q,T] x SI), then 

g(D(v N )) -± g(D(v N )) weakly inL p ([0,T] x ST). 

Now, we want to derive some bounds on the pressure. Let us start by the two 
dimensional case. 

By the Gagliardo Nirenberg inequality, we have 

,N\\ ^ /^ll„.JV||5 MV7„.JV|I5 



Thus, By Holder inequality, we get 

II JV-M4 < JV||2 \\X7v N \\ 2 

\\ V llL 4 ([0,T]xn) ° \\ V \\L°°(0,T,L 2 (n))\\ VV llL 2 ([0,T]xf7)- 

As p' < 2, then / <g> / is bounded in L p '([0, T] x ft). 
Since v N is divergence free, we get 

Ap N = divdiv [f{D{v N )) - v N <g> / + r w ) . 

The uniform boundedness of r w and f{D{v N )) in L p '([0,T] x ft) imply that p N is 
also bounded in L p '([0,T] x ft). Hence, we infer that 

p N -± p weakly inL p '([0,T] x Q). 

In the three dimensional case, the Gagliardo-Nirenberg inequality implies 

Since p > |, then Holder inequality leads to the fact that v N <S> v N is bounded in 
L p \[0,T] x Q). 
As above, we get 

p N ->> p weakly in L p ' ( [0, T] x SI) . 

Let us remark that at this stage, in the three dimensional case, we do not really need 
to have p > §. We need only to have p > 2 which lead to the fact that v N <g> v N is 
bounded in Lt(0, T, L p> (ST)) and thus ^ is bounded in L min ^''t)(0, T, L p '(Sl)). 

4.3. Strong convergence of v N in L 2 . We would like to prove that (v,T,p) is still 
a solution to (S). The difficulties appear when passing to the limit in the nonlinear 
terms 

v N .Vv N ; v N .Vr N ; f(D(v N )); g(D(v N )), r N .w(v N ) ~ w(v N ).r N . 
For the first two terms, we only need to prove that v N converges strongly to v in 
L 2 ([0, T] xSl). The bounds on the pressure and on the term v N ®v N imply that d t v N 
is bounded in L r (0, T, W~ 1,p ' (SI)), for some r e]l,+oo[. In addition, v N is bounded 
in L 2 (0, T, y 2 (Sl)). As V 2 (f2) is compactly embedded in L 2 (Sl) which is continuously 
embedded in W~ 1,p (ST), thanks to the Aubin-Lions lemma, up to a subsequence 

v N ^v strongly in L 2 ([0,T] x SI). 
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Hence, we get 

v N ® v N ->■ v ® v strongly in L x ([0,T] x O). 
Since v N is divergence free, then 

n n 

v ".Vt n = Tv?d,r N = Va(//) = div(vr) 



i=l i=l 

As converges strongly to u in L 2 ([0,T] x fi) and converges weakly to r in 
L 2 ([0,T] x Q), we deduce that 

v?t n -± ViT weakly mL l ([0, T] x O). 

Thus, it remains to prove that 



f(D(v N )) = f(D(v)), g{D{v N )) = g(D(v)) 
and 



T N .w(v N ) — w(v N ).T N = T.w(v) — w(v).T. 

If we prove that r N converges strongly to r in L 2 ([0,T] x Q), then our objective 
will be achieved. To get the strong convergence of r N , we will prove first that T N is 
equi-integrable in L 2 uniformly in t G [0,T]: 



lim sup sup / |T^| 2 Xri T iV| >M -i(i:E = 0. 
-^+°° n te[o,T]Jn 



4.4. Equi-integrability of t n in L 2 . We will show the equi-integrability of t n in 
L 2 uniformly in t G [0,T]. The idea is inspired from |28j . 

Let R > be fixed. We decompose t n into the sum of ip N and H N that solve 
respectively the following systems 

N id t tfj N + v N .Vip N + atfj N + i(j n .w(v n ) - w(v N ).ijj N = g(D(v N )), 
1 [ijj N (0,x) = t*(x)x { \ t n 1<r} , 



and 



N I d t H N + v N .VH N + aH N + H N .w(v N ) - w(v N ).H N 



Let M > 0, we have 



/ \T N \ 2 X{\r»\>M}dx = / \4> N + H N \ 2 X{\ T N\> M }dx 
Jn Jn 

< 2 / \4) N \ 2 X{\ T "\>M}dx + 2 / \H N \ 2 X{\ T N\> M} dx. 
Jn Jn 
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Energy estimates on H in the energy space L (Q) imply that 

/ \H N (t,x)\ 2 dx< I \H N (0,x)\ 2 dx=\\r N X{lT N l>R} \\l Hn) . 
Jn Jn 

Noticing that ^(0, .) G L P (Q), we should have a bound on ifj N in L°°(0,T, L P {Q)). 
We have 

l -d t \r? + a\r\ 2 + \v N .v\r\ 2 = = 

Multiplying this equation by p\ip N \ p ~ 2 , we get 

d t \ip N \ p + ap\^ N \ p + v N V\iP N \ p = p\ip N \ p - 2 g(D(v N )) : ^ N . 
As ft is bounded, integrating over Q, we deduce that 

d t [ \^ N (t,x)\ p dx <p [ ^{^x^-'lDiv^i^x^dx. 
Jn Jn 

Thanks to Holder inequality we get 

d t j \^ N (t,x)\ p dx< P (J \ip N {t, X )\ p dxy^ \D(v N )(t,x)\ p d X y. 

Multiplying this equation by ^ j Q \ip N (t, x)\ p dx^j v , we obtain 

d t U N (t,.)\\L, { n)< ( / \D(v N )(t,x)\ p dx)\ 

and thus by Holder inequality in time 

U N (t,.)\\ LP{n) < 11^(0, OII^+r^llD^llLPMxn). 
As D(v N ) is uniformly bounded in L p ([0,T] x Q), we get 

11^,011^(0) < ||^(0,.) + C. 
Applying Holder inequality, we obtain 

J \i> N \ 2 X{\r»\>M}dx < (f \ip N {t,x)\ p dxY [J X{\r^\>M}dx^ 

< C(\\ip N (0, .)\\ 2 LP (n) + C) (1 _2 J rAr |lL°°(oV,i 2 (o))- 

M p' 

Since t n is uniformly bounded in L°°(0, T, L 2 (fi)), then we obtain 

J \i> N \ 2 X{\r»\>M } dx < 67-^(11^(0, .)||L ( n) + C). 
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But, we have as p > 2 

11^(0, 011^(0) = [jr N (0,x)n {K]<R} dx 

< 



n 



T»(0,x)r 2 \T»(0,x)\% lToN{<R} dx 



< R p - 2 \\r N (0,.)\\l 2{n) . 
Since, we have \\P W^m-fh < 1> we deduce that 

11^(0, .)IU P( n) < if? 11^(0, .)||f 2(n) < ^||r ||| 2(n) . 
Finally, for a fixed R > 0, we get 

sup sup / \T N \ 2 X{\ T *\>M}dx < 2sup||r Af X{| r .jV|>fi}||i2ffi)+^ 1 2, (R J ^ A \\r \\ P L 2 {Q) +C). 
N te[0,T\JCl N M ( ?' 

As p > 2, we obtain 

(4.11) lim sup sup sup / \r N \ 2 x { \ T N\ >M} dx < 2 sup \\tq X{\ T » \>R}\\h(n)> 

m n te[o,T] Jn n 

which hold true for every i? g]0, +00 [. The fact that converges strongly to r in 
L 2 (Q) implies that (t^ 2 converges strongly in L 1 (fi) and thus weakly in L 1 (fi), and 
thus \t^\ 2 is equi-integrable in 

Consequently, let R —> +00 in inequality (14. lip , we infer that 

0<liminfsup sup / \r \ 2 X{\T N \>M}dx < limsupsup sup / \r N \ 2 X{\ T N \>M}dx < 0. 

This means that r N is equi-integrable in L 2 (Q) uniformly in t G [0,T]. 

4.5. Strong convergence of t n in L 2 . We will first focus on the system (S). 
As t n is equi-integrable in L 2 (f2) uniformly in t G [0, T], then |r — r| 2 converges 
weakly in Z^QO, T] x fi). So let 77 := \t n - r| 2 G L°°(0, T, L^O)). 
Notice that if 77 = 0, then 

It* - r| 2 -> strongly in #([0, T] x fi). 

Hence, our aim will be to show that 77 = 0. 
Multiplying by t n the equation satisfied by r , we get 

(4.12) (d t + v n V)\t n \ 2 + 2a\r N \ 2 = 2g{D{v N )) : r N . 

Notice here that the term v N .Vt n : t n has a sense since Vr N G L°°, v N G L 2 , t n G L 2 
for a fixe N G N. 

Let us introduce the unique a.e flow in the sense of DiPerna and Lions [10] of v , 
solution of 

d t X N {t, x) = v N (t, X{t, x)), X N (0, x) = x. 
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We also denote by X the a.e flow of v. 
Equation (I4.12p can be written as 

d t [\r N \ 2 (t,X N (t,x))} + 2a[\r N \ 2 (t,X N (t,x))] = 2^g(D(v N )) : r N )(t,X N (t,x)). 

Passing to the limit weakly, we get 

d t \T N \ 2 (t,X N (t,x)) + 2a\r N \ 2 (t,X N (t,x)) = 2{g{D{v N )) : r*) (t, X N (t, x)). 

The stability of the notion of a.e flow with respect to the weak limit of v N to v implies 
that X N (t,x) converges to X(t,x) in L\ oc and also that (X Ar (t)~ 1 )(a;) converges to 
{X(t)- l )(x) in L] oc . Thus, we get 

FFPW = \r»\*(t,X(t,x)) = (|r| 2 + V )(t, X(t, x)). 
Finally, we obtain 

(4.13) -(d t + v.V)(\t\ 2 + V ) + a(\r\ 2 + V ) = g{D{v N )) : r N . 

In addition, since v N .Vt n converges in D'([0, T] x Q) to i>.Vr, the equation satisfied 
by r is then 

d t r + ar + v.Vr + t n .w(v n ) -w(v N ).r N = g{D{v N )), 
which implies that 

(4.14) -(d t + v.V)\t\ 2 + a\r\ 2 + r N .w(v N ) - w(v N ).r N : r = g(D(v N )) : r. 

Subtracting (14.141) from (I4.13P gives the equation satisfied by r\ 

(4-15) 

d t g + 2ar] + div{r]v)-2r N .w{v N ) - w{v N ).r N : r = 2^g{D{v N )) : T N -g(D(v N )) : r). 

Notice that the term r/v makes sense in the sense of distributions in the two dimen- 
sional W 1,P (Q) C L°°(Q) and in the three dimensional case for p > 3, and also 
the term t n .w(v n ) — w(v n ).t n : r is not defined in the sense of distributions. To 
overcome this difficulty for any p, we use a renormalized form of (I4.15P by multiplying 
such equation by rjq^w, namely 

(d t + v.V)C + ^C = 77^(#W) : r N - g(D(v N )) : r 

+ t n .w(v n ) -w(v n ).t n : r) , 

where C = ^ e L°°([0,T] x Q). 
Let us remark that 

n := g(D(v N )) : - g(D(v»)) : r = (g(D(v N )) - g(D(v)j) : (t* - r) , 
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and 

T 2 := r N .w(v N ) -w(v n ).t n : r = (t n - t).(w(v n ) - w(v)) - (w(v N ) - w(v)).(t n - r) : r. 

As t n is equi-integrable in L 2 ([0,T] x Q) and w(v N ) is bounded in L 2 ([0,T] x fl), 
then (t n — t).(w(v n ) —w(v)) is equi-integrable in L 1 ([0, T] x Q) and thus converges 
weakly in L 1 ([0,T] x Q). Let z/ tiS be the Young measure generated by (t n , w(v N )), 
then we have 

(r N -r).(w(v N )-w(v)) = [ (a-T).(P-w(v))dv t , x (a,/3) 

J M nxn xM nxn 

By Cauchy-Schwartz inequality, we get 



(t"-t).(w(v")-w(v))\ < ( \a-r\ 2 du ttX (a,(3)Y 

^JM nxn xM nxn ' 

x ( / \(3-w(v)\ 2 dv t , x (a,(3)) 



-N ■ 



The equi-integrability of t n in L 2 ([0,T] x il) implies that 

V= |« - r| 2 rf^ i:E (a,/?). 

JM" x "xM" x " 

In addition, ^(v^)! 2 is equi-integrable in -^ 1 ([0, T] xfi) since it is bounded in La ([0, T] x 
f2) with p > 2, thus 

H^)-«^)l 2 = / \f3-w(v)\ 2 dv t:X (a,P). 

J M nxn xM nxn 

Consequently, we obtain 



\(t n - t).(w(v n ) - w(v))\ < ^q\w(v N ) - w(v)\ 2 \ 
and we remark that 



t n - t).(w(v n ) - w(v)) : r| < v - \t\\w{v n ) - w{v)\ 22 



(1 + V) 2,X '' v V ' v " • ' - (i + ?7 ) 

with G L°°([0, T] x ft), |r| G L 2 ([0, T] x O) and H^)-^)!^ G L 2 ([0, T] x 

ft). 

As r w is equi-integrable in L 2 ([0,T] x ft) and L>(u JV ) is bounded in L 2 ([0,T] x ft), 
then (t n - T).(g(D{v N )) - g(D(v))) is equi-integrable in ^([0,7] x ft). Let £ t):c be 
the Young measure generated by (r N , D(v N )), then we have 

Ti= f (a-T).(g(P)-g(D(v)))du t>x (a,P). 

J M nxn xM nxn 
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As the derivative of g is bounded, then g is Lipschitz and the same calculus as for 
the term T 2 imply that 

T^Cy/vlDiv^-Div)^. 

Finally, we get 

(dt + vV)( < C - ^ ^ (\D(v N ) -D(v)\ 2 * + \t\\w(v n ) - w(v)\ 2 *) . 

(1 + 77)^ v / 

Now, we should estimate \D(v N ) — D(v)\ 2 and \w(v N ) — w (v)\ 2 in terms of r\ with 
the help of the first equation in (S). More precisely, we will prove that 

(/(« " WW) '■■ (D(v) - D(v N )) = p - r) : (D(v N ) - D(v)). 
Multiplying the first equation in (S) by v N and using the identity 
div(A)u = div(Au) - A : D(u), A e M nxn , u e R n , 

we get 
(4.16) 

l -d t \v N \ 2 -div{f{D{v N ))v N )+f{D{v N )) : D{v N )+v N .Vv N .v N +div(p N v N ) = div{r N v N )-T N : D(v N ). 

The strong convergence of v N to v in L 2 ([0,T] x Q) implies that \v N \ 2 converges 
strongly in L l ([Q,T] x 9) to \v\ 2 . 

In addition, thanks to the weak convergence of t n to r in L 2 ([0,T] x Q), we can 
conclude that 

r N v N ->> Tf weakly in £*([(), T] x ft). 
But, v N is uniformly bounded in L p (0,T, V p ) and is uniformly bounded in 

L p '(0, T, V4 7 ~ 1 ' p '(fi)) in the two dimensional case and L min fr'>i)(0, T, W' 1 ^'^)) in the 
three dimensional case. Since V p is compactly embedded into L p which is continuously 
embedded into W~ 1,p ' (Q), Rellich-Kondrachov Theorem implies that 

v N -> v strongly in L p ([0, T] x ft). 

As f(D(v n )) is uniformly bounded in L p '([0,T] x ft), we obtain 

f(D(v n ))v N -± f(D(v n ))v weakly inL^CT] x ft). 

In the two dimensional case, p N and v N ®v N are uniformly bounded in L p '([0, T] x ft). 
Thus, one gets 

p N v N -± pv weakly in L 1 ([0, T] x ft). 
Notice that thanks to the divergence free condition on v N , 

v N .Vv N .v N = l -div{\v N \ 2 v N ), 

and thus 

\v N \ 2 v N -± \v N \ 2 v = \v\ 2 v weakly inL^O,! 1 ] x ft). 
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In the three dimensional case, for p > |, p N and v N <S> v N are uniformly bounded in 
U ([0, T] x ft). Notice that at this stage appears the bound on p when the dimension 
is three. In fact, we have 



p' 3p' 

< Ih^ll 7 \\Vv N \\~ 

^ llL 2 (n)H vu \\LP(n)- 



Thus, it suffices to have ^- < p that is p > 
Taking the weak limit of f !4.16p . we obtain 
(4.17) 
1 



^d t \v\ 2 -div(f(D(v N ))v)+f(D(v N )) : D(v N )+v.Vv.v+div(pv) = div(rv)-T N : D(v N ) 

On the other hand, we take the weak limit of the first equation in (S) knowing that 
v N .Vv N converges to v.Vv 



d t v + v.Vv - div(f(D(v N ))) +Vp = divr, 

which is multiplied by v leads to 
(4.18) 



1- ,2 



d t \v\ 2 -div(f(D(v N ))v) + f(D(v N )) : D(v)+v.Vv.v++div(pv) = div(Tv)-r : D(v). 
Subtracting (I4.18P from f)4.17p . we get 



f(D(v N )) : D(v N ) - f(D(v N )) : D(v) = -r N : D(v N ) + r : D(v). 
Or equivalently 



(/(£>(«")) - f(D(v))) : (D(v) - £>(„")) = (r N - r) : (£>(„") - £>(«)). 

For the system (S2), under the hypotheses on / in Theorem I3.6[ since / is monotone, 
we deduce that 



(t n - r) : (D(v N ) - D(v)) < 0, 
and thus, the equation satisfied by £ is reduced to 

2 

which implies that ( = rj = 0. 

Now, for the system (Si), under the strong monotonicity condition on / in Theorem 
13. 5\ we deduce that 



(d t + vV)( < rr-^^ - r) : (D(v N ) - D(v)) < 0, C(0, .) = 0, 



\D(v N ) -D(v)\ 2 < \(t n -t) : (D(v N ) - D(v))\ < ^/rj\D(v N ) - D(v) 



2-> 



< -r,+ -\D(v^-D(v)\ 2 . 
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Thus, we get 



\D(v N ) - D(v)\ 2 < rj, 

and 

(d t + v.V)(<C(, 

which implies thanks to Gronwall's lemma and to the fact that £(0, .) = that 

C = v = o. 

Notice that if the quadratic term t n .w(v ) — w(v N ).r N is present, then the difficulty 
is how to estimate \w(v N ) — w(v)\ 2 in terms of rj. 

4.6. Weak limits of non linear terms involving D(v N ). Under the hypotheses 
on / in Theorem 13 .5} we have 



\D(v N ) - D(v)\ 2 < r, = 0. 

As \D(v N ) — D(v)\ 2 is bounded in ([0, T] x Q) with p > 2, then it is equi-integrable 
in ^([OjT] x Q), and hence 

D(v N ) -> D(v) strongly inL 2 ([0,T] x Q). 

Vitali's lemma implies that 



f(D(v N )) = f(D(v)), g(D(v»)) = g(D(v)). 



Let us now focus on the system (S2). Under the hypotheses on / in Theorem 13 .6[ we 
have two ways to get f(D(v N )). 

Let us begin by the simplest way. Let /i t;X be the Young measure generated by the 
sequence D(v N ) and let G(X,t,x) : = (/(A) - f(D(v)(t,x))) : (A - D(v)(t,x)) > 0. 
Remark that 

I G(\,t, x)p, tiX (X) — [ lim fiif ^ * , v dfx tjX (A) . 

J M nxn J M nxn 1 ~T OLr^A, t, 

By Fatou's lemma, we get 



G(\,t,x)dfAt,x{ty < liminf / — 

<5— >0 J MnXn 1 



G{X,t,x) 
+ 5G(X,t,x) 



d/j, t ,x(X). 



Notice that 



G(D(v N ),t,x) ^1 G(D(v N ),t,x) ^1 



l + - 51 + G(D(v N ),t, x) ~ <f 

Hence, for a fixed 5 e]0, 1], the sequence j^^fe^i^) is bounded in L°°([0,T] x Q) 
and thus equi-integrable in L 1 ([0,T] x Q). 
We can claim that 



G(D(v»),t, X ) _ f G{\,t,x) ^ (A)>vfe]0il] . 



1 + 5G(D(v N ),t, x) J M nxn 1 + 6G(X,t,x 
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Consequently, we get 



/ 



G(D(v N ),t,x) 



But, thanks to the strong convergence of t n to r in L 2 ([0,T] x fi) and to the weak 
convergence of D(v N ) to U('y) in L 2 ([0,T] x we deduce that 



G(X,t,x)d^ x (X) < G(D(v N ),t,x) = -(r* - r)(D(v") - D(v)) = 0. 

M nxn 

This means that the div-curl inequality used in [8] and [9] , which is the key ingredient 
to handle the limits of non-linear terms, remains true in our case and enables us to 
get /(£>K)) = f(D(v)). 

The second way is to proceed as in [S] and [H] to get the div-curl inequality. 

As / is monotone, then the negative part of (f(D(v N )) — f(D(v))) : (D(v N ) — D(v)) 

is null and thus weakly relatively compact in L 1 . Applying Theorem 13 A\ we infer 

that 

(/(A) - /(£>(«))) : (A - D(v))dfi SyX (X)dxds 




o J n J m 



r nXn 



< liminf / / (f(D(v N )) - f(D(v))) : (D(v N ) - D(v))dxds 

rw+oo J Q J n 



< limsup / / (f(D(v N )) - f(D(v))) : (D(v N ) - D{v))dxds 

n->+oo Jo Jfl 



n— >+oo 



t r rt 

N 



< -liminf / / f(D(v)) : D(v N )dxds + / / f(D(v)):D(v)dxds 
'o in Jo Jo. 



+ limsup/" [ f(D(v N )) : D(v N )dxds 
- liminf / / f(D(v N )) : D(v)dxds. 



The weak convergence in L p of D(v N ) to -D(w) and the fact that f(D(v)) belongs to 
L p imply that 



n— >+oo 



t r rt 

N 



liminf / / f(D(v)) : D(v ls )dxds = / / : D(v)dxds. 

'o Jn Jo Jn 



Moreover, we have 



ft r rt r 

liminf / / f(D(v N )) : D(v)dxds = f(D(v N )) : D(v)dxds. 

'o Ju Jo Jn 



n— >+oo 



GLOBAL WEAK SOLUTIONS ... 23 

N 



Testing the first equation in (S) by v and integrating over [0, t] x f2, we get 

f [ f(D(v N )) : D(v N )dxds =11 (--d t \v N \ 2 - r N : D{v N ))dxds 
Jo Jq Jo Jq 2 

< ~J*I tN -- D(v N )dxds+ \\\v»\\h m - \\\v N (t,.)\\h { ny 

The same arguments as for the classical Navier-Stokes equation lead to 

v N (t,.) ->> v(t,.) weakly in L 2 ,W G [0,T]. 
Thanks to the lower semi-continuity of the norm, we have 

\\v(t, -)IU 2 (n) < liminf \\v N (t, -)IU 2 (n)- 



Using the fact that Vq converges strongly to v in L , we obtain 

limsup f [ f(D(v N )) : D(v N )dxds < l\\v \\ 2 L2{Q) + limsup / [ -t n : D(v N )dxds 
n-^+oo Jo Jq 2 n^+oo Jo Jq 



Knowing that 



-d t \v\ 2 -div(f(D(v N ))v) + f(D(v N )) : D(v) +v.Vv.v + div(pv) =div(rv)-T : D(v), 
we get by integration over [0,t] x Q 

f(D(v»)) : D{v)dxds = ~\\v(t, .)\\ 2 LHn) + hv \\ 2 L2{n) - f f r : D(v)dxds. 




J fi 

Finally, we get 

1 [ [ (/(A) - /(£>(«))) : (A - D(v))d^ x (X)dxds 
o Jq J M nxn 

< 



liminf / / (f(D(v N )) - f(D(y))) : (D(v N ) - D(v))dxds 
< limsup / / -r w : D(v N )dxds + / / r : D(v)dxds. 

Since r N converges strongly to r in L 2 and D(v N ) converges weakly in L 2 to D(v), 
then r N : D(v N ) converges weakly in L 1 to r : D{v), particularly 

lim f f t n : D(y N )dxds = limsup / f t n : D(v N )dxds = [ It: D(v)dxds. 
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Hence 

(f(X) - f(D(v))) : (X-D(v))d^ x (X)dxds 




JCiJM nXn 



< liminf / / (f(D(v N )) - f(D(v))) : (D(v N ) - D(v))dxds 



< o. 

This ends the proofs of Theorems l3.5H3.6i 



Acknowledgments. We are very grateful to Professor Nader Masmoudi for 
interesting discussions around the questions dealt with in this paper. 

References 

[1] H. W. Alt, Lineare Funktionalanalysis (2nd edition), Springer- Verlag, Berlin- Heidelberg- New 
York (1992). 

[2] C. Andrade, Viscosity of liquids, Nature, 125, 309-310, (1930). 

[3] N. Arada, A. Sequeira, Strong steady solutions for a generalized Oldroyd-B model with shear- 
dependent viscosity in a bounded domain, Math. Mod. Meth. Appl. S, 13(9), 1-21, (2003). 

[4] H. Brezis, Analyse fonctionnelle Theorie et applications, Masson, Paris, (1983). 

[5] P. W. Bridgman, The physics of high pressure, MacMillan, New York, (1931). 

[6] M. Bulicek, M. Majdoub, J. Malek, Unsteady flows of fluids with pressure dependent 
viscosity in unbounded domains, Nonlinear. Anal. Theor, 11(5), 3968-3983, (2010). 

[7] J. Y. Chemin, N. Masmoudi, About lifespan of regular solutions of equations related to 
viscoelastic fluids, SIAM J. Math. Anal, 33(1), 84-112 (electronic), (2001). 

[8] P. Dreyfuss, N. Hungerbuhler, Results on a Navier-Stokes system with applications to 
electrorheological fluid flow, Int. J. Pure. Appl. Math, 14(2), 241-271, (2004). 

[9] P. Dreyfuss, N. Hungerbuhler, Navier-Stokes systems with quasimonotone viscosity tensor, 
Int. J. Differ. Equ. Appl, 9(1), 59-79, (2004). 
[10] R. J. DiPerna, P. L. Lions, Ordinary differential equations, transport theory and Sobolev 

spaces, Invent. Math, 98(3), 511-547, (1989). 
[11] E. Fernandez- Cara, F. Guillen, R. R. Ortega, Some theoretical results for viscoplastic 

and dilatant fluids with variable density, Nonlinear Anal, 28(6), 1079-1100, (1997). 
[12] E. Fernandez- Cara, F. Guillen, R. R. Ortega, Some theoretical results concerning non- 
Newtonian fluids of the Oldroyd kind, Ann. Scuola Norm. Sup. Pisa CI. Sci (4), 26(1), 1-29, 
(1998). 

[13] E. Fernandez- Cara, F. Guillen, R. R. Ortega, The mathematical analysis of viscoelastic 

fluids of the Oldroyd kind, (2000). 
[14] C. Guillope, J. C. Saut, Existence results for the flow of viscoelastic fluids with a diffferential 

constitutive law, Nonlinear Anal, 15(9), 849-869, (1990). 
[15] C. Guillope, J. C. Saut, Global existence and one- dimensional nonlinear stability of shearing 

motions of viscoelastic fluids of Oldroyd type, RAIRO Model. Math. Anal. Numer, 24(3), 369- 

401, (1990). 

[16] R. R. Huilgol, Continuum mechanics of viscoelastic liquids, Hindusthan Publishing Corpora- 
tion, Delhi, (1975). 



GLOBAL WEAK SOLUTIONS 



25 



[17] N. Hungerbuhler, A refinement of Ball's Theorem on Young measures, New York J. math, 
3, 48-53, (1997). 

[18] R. Kupferman, C. Mangoubi, E. S. Titi, A Beale-Kato-Majda breakdown criterion for an 
Oldroyd-B fluid in the creeping flow regime, Math. Sci, 6(1), 235-256, (2008). 

[19] Z. Lei, Rotation- Strain Decomposition for the Incompressible Viscoelasticity in Two Dimen- 
sions, http : //arxiv . org/ abs/1204 . 5763 . 

[20] Z. Lei, Global existence of classical solutions for some Oldroyd-B model via the incompressible 
limit, Chinese Ann. Math. Ser. B, 27(5), 565-580, (2006). 

[21] Z. Lei, On 2D viscoelasticity with small strain, Arch. Ration. Mech. Anal, 198(1), 13-37, 
(2010). 

[22] Z. Lei, C. Liu, Y. Zhou, Global existence for a 2D incompressible viscoelastic model with 

small strain, Commun. Math. Sci, 5(3), 595-616, (2007). 
[23] Z. Lei, C. Liu, Y. Zhou, Global solutions for incompressible viscoelastic fluids, Arch. Ration. 

Mech. Anal, 188(3), 371-398, (2008). 
[24] Z. Lei, N. Masmoudi, Y. Zhou, Remarks on the blowup criteria for OLdroyd models, J. 

Differ. Equations, 248(2), 328-341 , (2010). 
[25] Z. Lei, Y. Wang, Global solutions for micro-macro models of polymeric fluids, J. Differ. 

Equations, 250(10), 3813-3830, (2011). 
[26] Z. Lei, Y. Zhou, Global existence of classical solutions for the two-dimensional Oldroyd model 

via the incompressible limit, SIAM J. Math. Anal, 37(3), 797-814 (electronic), (2005). 
[27] J. L. Lions, Quelques methodes de resolution des problemes aux limites non lineaires, Dunod, 

Paris, (1969). 

[28] P. L. Lions, N. Masmoudi, Global solutions for some Oldroyd models of non-Newtonian 
flows, Chinese Ann. Math. Ser. B, 21(2), 131-146, (2000). 

[29] J. Malek, J. Necas, K. R. Rajagopal, Global analysis of the flows of fluids with pressure- 
dependent viscosities, Arch. Ration. Mech. Anal, 165 (3), 243-269, (2002). 

[30] J. Malek, J. Necas, M. Rokyta, M. Ruzicka, Weak and Measure-Valued Solutions to 
Evolutionary PDEs, Chapman-Hall, (1996). 

[31] J. Malek, K. R. Rajagopal, Mathematical properties of the solutions to the equations gover- 
ing the flow of fluid with pressure and shear rate dependent viscosities. In: Handbook of Mathe- 
matical Fluid Dynamics, Handb. Differ. Equ. Elsevier /North- Holland, Amsterdam 4, 407-444, 
(2007). 

[32] F. Murat, J. M. Ball, Remarks on Chacon's biting lemma, Proc. Amer. Math. Soc, 107(3), 
655-663, (1989). 

[33] T. Rubicek, A generalization of the Lions-Temam compact imbedding theorem, Casopis pest. 

mat, 115, 338-342, (1990). 
[34] W. R. Schowalter, Mechanics of non- Newtonian fluids, Pergamon Press, Oxford, (1978). 
[35] A. Z. Szeri, Fluid Film Lubrication: Theory and Design, Cambridge University Press, (1998). 
[36] G. B. Thurston, Viscoelasticity of human blood, Biophys. J, 12, 1205-1217, (1972). 
[37] W. Walter, Differential and integral inequalities, Springer- Verlag, Berlin-Heidelberg-New 

York, (1970). 

[38] K. Yosida, Functional analysis, Springer- Verlag, Berlin-Gottingen- Heidelberg, (1965). 
[39] E. Zeidler, Nonlinear Functional analysis II/A-Linear Monotone Operators, Springer- Verlag, 
Berlin-Heidelberg- New York, (1990). 



26 



OLFA BJAOUI AND MOHAMED MAJDOUB 



University of Tunis ElManar, Faculty of Sciences of Tunis, Department of Math- 
ematics, ElManar 2092, Tunisia 

E-mail address: bjaouiolfa@yahoo.fr 

University of Tunis ElManar, Faculty of Sciences of Tunis, Department of Math- 
ematics, ElManar 2092, Tunisia 

E-mail address: mohamed .majdoub@f st .rmi.tn 



